Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Mouse cellular cementum is highly dependent on growth hormone status

Cementum is known to be growth-hormone (GH)-responsive, but to what extent is unclear. This study examines the effects of extremes of GH status on cementogenesis in three lines of genetically modified mice; GH excess (giant), GH antagonist excess (dwarf), and GH receptor-deleted (GHR-KO) (dwarf). Age-matched mandibular molar tissues were processed for light microscope histology. Digital images of sections of first molar teeth were captured for morphometric analysis of lingual root cementum. Cross-sectional area of the cellular cementum was a sensitive guide to GH status, being reduced nearly 10-fold in GHR-KO mice, three-fold in GH antagonist mice, and increased almost two-fold in giant mice (p

Semliki Forest virus and Kunjin virus RNA replicons elicit comparable cellular immunity but distinct humoral immunity

RNA replicons offer a number of qualities which make them attractive as vaccination vectors. Both alphavirus and flavivirus replicon vaccines have been investigated in preclinical models yet there has been little direct comparison of the two vector systems. To determine whether differences in the biology of the two vectors influence immunogenicity, we compared two prototypic replicon vectors based on Semliki Forest virus (SFV) (alphavirus) and Kunjin virus (KUN) (flavivirus). Both vectors when delivered as naked RNAs elicited comparable CD8+ T cell responses but the SFV vectors elicited greater humoral responses to an encoded cytoplasmic antigen beta-galactosidase. Studies in MHC class II-deficient mice revealed that neither vector could overcome the dependence of CD4+ T cell help in the development of humoral and cellular responses following immunization. These studies indicate that the distinct biology of the two replicon systems may differentially impact the adaptive immune response and this may need to be considered when designing vaccination strategies. (c) 2005 Elsevier Ltd. All rights reserved.

Understanding alternative splicing: Towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

Assessing the cellular transmembrane electrical potential difference on the hepatic uptake of palmitate

Understanding the driving forces for the hepatic uptake of endogenous and exogenous substrates in isolated cells and organs is fundamental to describing the underlying hepatic physiology/pharmacology. In this study we investigated whether uptake of plasma protein-bound [H-3]-palmitate across the hepatocyte wall is governed by the transmembrane electrical potential difference (PD). Uptake was studied in isolated hepatocytes and isolated perfused rat livers (IPL). Protein-binding and vasoactive properties of the different perfusates were determined using in vitro heptane/buffer partitioning studies and the multiple indicator dilution (MID) technique in the IPL, respectively. Altering hepatocyte PD by perfusate ion substitution resulted in either a substantial depolarization (-14 +/- 1 mV, n = 12, mean +/- S.E., substituting choline for Na+) or hyperpolarization (-46 +/- 3 mV, n = 12, mean +/- S.E., substituting nitrate for Cl-). Perfusate ion substitution also affected the equilibrium binding constant for the palmitate-albumin complex. IPL studies suggested that, other than with gluconate buffer, hepatic [H-3]-palmitate extraction was not affected by the buffer used, implying PD was not a determinant of extraction. [H-3]-Palmitate extraction was much lower (p < 0.05) when gluconate was substituted for Cl- ion. This work contrasts with that for the extraction of [H-3]-alanine where hepatic extraction fraction was significantly reduced during depolarization. Changing the albumin concentration did not affect hepatocyte PD, and [H-3]-palmitate clearance into isolated hepatocytes was not affected by the buffers used. MID studies with vascular and extravascular references revealed that, with the gluconate substituted buffer, the extravascular volume possibly increased the diffusional path length thus explaining reduced [H-3]-palmitate extraction fraction in the IPL.

Are insect immune suppressors driving cellular uptake reactions?

Many insect parasitoids that deposit their eggs inside immature stages of other insect species inactivate the cellular host defence to protect the growing embryo from encapsulation. Suppression of encapsulation by polydnavirus-encoded immune-suppressors correlates with specific alterations in hemocytes, mainly cytoskeletal rearrangements and actin-cytoskeleton breakdown. We have previously shown that the Cotesia rubecula polydnavirus gene product CrV1 causes immune suppression when injected into the host hemocoel. CrV1 is taken up by hemocytes although no receptors have been found to bind the protein. Instead CrV1 uptake depends on dimer formation, which is required for interacting with lipophorin, suggesting a CrV1-lipophorin complex internalisation by hemocytes. Since treatment of hemocytes with oligomeric lectins and cytochalasin D can mimic the effects of CrV1, we propose that some dimeric and oligomeric adhesion molecules are able to cross-link receptors on the cell surface and depolymerise actin by leverage-mediated clearance reactions in the hemolymph.

Kidney side population reveals multilineage potential and renal functional capacity but also cellular heterogeneity

Side population (SP) cells in the adult kidney are proposed to represent a progenitor population. However, the size, origin, phenotype, and potential of the kidney SP has been controversial. In this study, the SP fraction of embryonic and adult kidneys represented 0.1 to 0.2% of the total viable cell population. The immunophenotype and the expression profile of kidney SP cells was distinct from that of bone marrow SP cells, suggesting that they are a resident nonhematopoietic cell population. Affymetrix expression profiling implicated a role for Notch signaling in kidney SP cells and was used to identify markers of kidney SP. Localization by in situ hybridization confirmed a primarily proximal tubule location, supporting the existence of a tubular niche, but also revealed considerable heterogeneity, including the presence of renal macrophages. Adult kidney SP cells demonstrated multilineage differentiation in vitro, whereas microinjection into mouse metanephroi showed that SP cells had a 3.5- to 13-fold greater potential to contribute to developing kidney than non-SP main population cells. However, although reintroduction of SP cells into an Adriamycin-nephropathy model reduced albuminuria:creatinine ratios, this was without significant tubular integration, suggesting a humoral role for SP cells in renal repair. The heterogeneity of the renal SP highlights the need for further fractionation to distinguish the cellular subpopulations that are responsible for the observed multilineage capacity and transdifferentiative and humoral activities.

Focusing in on structural genomics: The University of Queensland structural biology pipeline

The flood of new genomic sequence information together with technological innovations in protein structure determination have led to worldwide structural genomics (SG) initiatives. The goals of SG initiatives are to accelerate the process of protein structure determination, to fill in protein fold space and to provide information about the function of uncharacterized proteins. In the long-term, these outcomes are likely to impact on medical biotechnology and drug discovery, leading to a better understanding of disease as well as the development of new therapeutics. Here we describe the high throughput pipeline established at the University of Queensland in Australia. In this focused pipeline, the targets for structure determination are proteins that are expressed in mouse macrophage cells and that are inferred to have a role in innate immunity. The aim is to characterize the molecular structure and the biochemical and cellular function of these targets by using a parallel processing pipeline. The pipeline is designed to work with tens to hundreds of target gene products and comprises target selection, cloning, expression, purification, crystallization and structure determination. The structures from this pipeline will provide insights into the function of previously uncharacterized macrophage proteins and could lead to the validation of new drug targets for chronic obstructive pulmonary disease and arthritis. (c) 2006 Elsevier B.V. All rights reserved.

Parallel implementation of stochastic simulation for large scale cellular processes

Experimental and theoretical studies have shown the importance of stochastic processes in genetic regulatory networks and cellular processes. Cellular networks and genetic circuits often involve small numbers of key proteins such as transcriptional factors and signaling proteins. In recent years stochastic models have been used successfully for studying noise in biological pathways, and stochastic modelling of biological systems has become a very important research field in computational biology. One of the challenge problems in this field is the reduction of the huge computing time in stochastic simulations. Based on the system of the mitogen-activated protein kinase cascade that is activated by epidermal growth factor, this work give a parallel implementation by using OpenMP and parallelism across the simulation. Special attention is paid to the independence of the generated random numbers in parallel computing, that is a key criterion for the success of stochastic simulations. Numerical results indicate that parallel computers can be used as an efficient tool for simulating the dynamics of large-scale genetic regulatory networks and cellular processes